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ABSTRACT 

Context. Turbulence is thought to be a key driver of the evolution of protoplanetary disks, regulating the mass accretion process, the 
transport of angular momentum, and the growth of dust particles. 

Aims. We intend to determine the magnitude of the turbulent motions in the outer parts (> 100 AU) of the disk surrounding DM Tau. 
Methods. Turbulent motions can be constrained by measuring the nonthermal broadening of line emission from heavy molecules. We 
used the IRAM Plateau de Bure interferometer to study emission from the CS molecule in the disk of DM Tau. High spatial (1.4 x 1'.') 
and spectral resolution (0.126 km s ) CS J=3-2 images provide constraints on the molecule distribution and velocity structure of the 
disk. A low sensitivity CS 1=5-4 image was used in conjunction to evaluate the excitation conditions. We analyzed the data in terms 
of two parametric disk models, and compared the results with detailed time-dependent chemical simulations. 

Results. The CS data confirm the relatively low temperature suggested by observations of other simple molecules. The intrinsic 
linewidth derived from the CS 1=3-2 data is much larger than expected from pure thermal broadening. The magnitude of the derived 
nonthermal component depends only weakly on assumptions about the location of the CS molecules with respect to the disk plane. 
Our results indicate turbulence with a Mach number around 0.4 - 0.5 in the molecular layer. Geometrical constraints suggest that this 
layer is located near one scale height, in reasonable agreement with chemical model predictions. 

Key words. Stars: circumstellar matter - planetary systems: protoplanetary disks - individual: DM Tau - Radio-lines: stars 



1. Introduction 



Turbulence is thought to be a major actor throughout the life- 
time of protoplanetary disks, by regulating the mass accretion 
and angular momentum transport, planetesimal formation and 

proces ses relevant for gas and grain 
The magneto-rotational in- 



planet migration, an d mixing proces 
chemistry (see, e.g. lHennindl2008l) 



stability (MRI, iBalbus & Hawlevlll99ll) is currently the most 
widely accepted theory for an origin of the effective turbu- 
lent viscosity. The ionization fraction in the outer disk is suf- 
ficient to allow for MRI to operate, but there may be a "dead" 
laminar zone in the inner (< 10 AU) disk inte rior, efficiently 
shielded from ionizing radiation dGammidl 19961) . Its extent de- 
pends on a variety of factors including magnetic field morphol- 
ogy, X-ray activity, ab undance of metals, and d ust properties 
Ilgner & Nelsonll2006l). The 2D disk models by iKeller & Gaill 



2004 and iTscharnuter & Gaill J2007h show the possible pres- 



ence of large-scale circulation (lUrpinll 19841: iRegev & Gitelmanl 
|2002|) . but such motions would be suppressed by turbulence on 
the scale of the disk thickness. These meridional flows may 
be artifacts of the 2D-hydrodynamical disk models with fixed 
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a viscosity. More realistic 3D MHD simulations do not show 
preferential direc tions of the large-scale gas-flows (see e.g., 
iFlock et al.ll20TTh 

Turbulence may also provide a source of relative ve- 
lociti es for grain coagulati on to occur (e.g. IVoelk et all 
119801: iBeckwith etaf] 120001) . The interplay between gas 
and solid particle motion, would also lead to high local 
"dust" overdensi ties and planetesimal forma tion by gravita- 
tional collapse dJohansen et al.l 120071 1201 ll) . The turbulent 
state may have a strong imp act on the migration efficiency 
dPapaloizou & Terqueml l2006t lOishi et al.1 120071: lUribe et al.l 
1201 ll) . Finally, models of gas-phase chemistry including sur- 
face chemisny on grains are sensitive to the level of tur- 
bulent mixing (e.g. Illgner et al.l 12004 lllgner & Nelson! 120061: 
ISemenov & Wiebell201 H) w hich may partially explain the pres- 
ence of cold CO in disks dDartois et al.l 120031: ISemenov et al.l 
l2006HAikawall2007l:lHersant et al.ll2009l 

Despite this enormous importance of turbulence for disks, 
even its magnitude is poorly constrained. Most constraints on 
the anomalous viscosity parameter a come from measured ac- 
cretion rates on to the star, which is mos tly determined by inner 
disk properties. iHueso & Guillotl (120051) have attempted to con- 
strain the overall value of the a parameter by simulations of the 
viscous evolution of disks as a function of ti me. More recent 
results on viscous spreading of dust disks from lGuilloteau et al.l 
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(1201 lh suggest a decreases with time, with values less than 0.001 
near 100 AU after a few Myr. 

Direct measurements of the nonthermal gas moti ons result- 
ing fro m mesoscopic scale turbulence remain rare. ICarr et ail 
d2004l) report supersonic turbulence from CO overtone bands, 
but lower linewidths for HwO in S VS 13, suggesting strong turbu- 
lence in the upper layers of the inner disk . On th e other hand, for 
the outer disk of DM Tau, iDartois et alj d2003l) find marginally 
larger local velocity dispersion in 13 CO than in 12 CO, suggest- 
ing lower levels of turbulence in the upper layers than in the disk 
plane. The derived broadening, 0.07 to 0.14 km.s -1 , is subsonic, 
but its exact value is difficult to assess because of the contribu- 
tion of the thermal component. iPietu et alJ (12007) find similar 
results for MWC 480 and LkCa 15. More recenflv. lHughes et alj 
(2011) have obtained a low value for TW Hya (< 0.04 km.s '), 
but up to 0.3 km.s" 1 for the Herbig Ae star HD 163296 from CO 
J=3-2 observations. 

However, although 12 CO is an interesting indicator, because 
of the significant opacities of the lower lying transitions, it only 
probes a thin upper layer in the disk, where strong temperature 
gradients may affect the apparent linewidth. To build a com- 
plete picture of turbulence in disks, it is thus essential to con- 
strain the turbulence from other (optically thinner) tracers. CO 
isotopologues and other molecules like HCO + and CN yield 
typical linewidths around 0.18 km.s -1 at 300 AU, but attribut- 
ing this to thermal or nonthermal motions re mains difficult, as 
their thermal width at 30 K is 0.14 km.s -1 dPietu et al J 120071: 
IChanillon et afll2012bh . Moreover, current chemical models fail 
to explain their vertical location in the molecular layer since CN, 
CCH o r HCN appear to be in a colder molecular la yer than pre- 
dicted dHenning et alJl2010HChapillon et al"]l2012bl) . 

It is then important to reduce the uncertainty (and possible 
bias) introduced by the thermal contribution to the linewidth. 
Among molecules detected in disks, only CS is both abundant 
enough and sufficiently heavier than the previously cited ones. 

We report here on the study of the CS molecule in the DM 
Tau disk. Observations are presented in Section [2 disk proper- 
ties are derived in Section[3j and we discuss the results in Section 

El 



2. Observations 

Observations were carried out with the IRAM Plateau de Bure 
interferometer. We observed two different transitions of CS, the 
J=3-2 line at 146.969047 GHz and the J=5-4 line at 244.935609 
GHz. 

For the J=3-2 transition, we used the BC configuration 
(project RB8F). The B configuration, with baselines up to 408 m 
was observed on 18 Mar 2008, while the more compact C config- 
uration was observed on 15 Nov 2008. Phase noise ranged from 
15° up to 55° on the longest baselines. Amplitude calibration 
was done using nearby quasars and led to an rms level of about 
5%. Absolute flux calibration was done using MWC 349 as a ref- 
erence. The angular resolution is 1 .4" x 1 .0". The correlator pro- 
vides only a limited spectral resolution with a channel spacing 
of 0.08 km.s -1 , but with an effective resolution of 1.59 times the 
channel spacing (i.e. 0.126 km.s -1 ) owing to Welch apodization 
of the Fourier spectrum. At this resolution, the brightness sensi- 
tivity is 8.8 mJy/beam, which corresponds to 0.35 K at this angu- 
lar resolution. Channel maps of CS(3-2) obtained after smooth- 
ing in velocity to 0.16 km.s -1 are presented in FigQ] but our 
analysis uses the unsmoothed data. 

The CS(5-4) data were obtained on 11 Nov 2004, 17 Nov 
2004, and 03 Apr 2005 (Project OA45). The natural weighting 



CS (5-4) 



Fig. 2. Spectrum of CS J=5-4 line towards DM Tau. The angular 
resolution is 3.5" and the spectral resolution 0.25 kms -1 . The 
curve is a best-fit Gaussian to emphasize the detection. 



spatial resolution is 2.6" x 1.8" at PA 85°. We smoothed the 
data to a spectral resolution of 0.25 km.s -1 for the analysis. The 
brightness sensitivity is 90 mJy/beam, corresponding to 0.4 K. 
The CS(5-4) is barely detected (at the 5<x level): a spectrum to- 
wards DM Tau obtained at 3.5" resolution is shown in Figf2] 
The data is com patible with the int egrated spectrum obtained 
with the 30-m bv lDutrev et al.l d 19971) and a source size of about 
5", with an integrated line flux 0.6 ± 0.1 Jy.km.s -1 . 

The relative calibration accuracy between the two frequen- 
cies is better than 10 %. It is based on the precise spectral index 
of 0.6 for the reference flux calibration source MWC 349. 



3. Results 

3. 1 . Principle of the method 

Line formation in a Kepler ian disk is strongly con strained by the 
velocity gradient (see, e.g. lHorne & Marshll 19861 in the context 
of cataclysmic binaries). The line of sight velocity (in the system 
rest frame) is 



Vo bs (r,0)= jGMJr sin i cos 6 



(1) 



where r, 9 are the cylindrical coordinates in the disk plane. The 
locii of isovelocity are given by 



r{0) = (GM,/7 2 bs ) sin 2 ; cos 2 6. 



(2) 



With a finite local linewidth Av (assuming rectangular line shape 
for simplification), the line at a given velocity V b s originates in 
a region included between r,{9) and r s {0) : 



n{0) 

r s {6) — min 



(V obs + Av/2) 2 



sin 2 i cos 2 9 



GM„ 



(V obs - Av/2) 2 



sin 2 i cos 2 9 



(3) 
(4) 



Figure [3] indicates the regions of equal projected velocities for 
six different values: V b s > \'d, V bs = Vj, V bs < Vd, and their 
symmetric counterpart at negative velocities, where Vd 



Vd 



y/GM t /R out sin/ 



(5) 



is the projected velocity at the outer disk radius R out . Figure[3]and 
Eq|5]thus show that the basic shape of the emission as a func- 
tion of velocity depends only on inclination i, while the stellar 
mass gives the scaling factor for the velocity scale. The intensity 
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Right Ascension 

Fig. 1. Channel maps of the CS J=3-2 emission towards DM Tau. The angular resolution is 1.4" x 1.0" and the spectral resolu- 
tion 0.16 kms~'. Contour spacing is 10 mJy/beam, corresponding to 2<x and 0.4 K brightness. The cross indicates the position, 
orientation, and aspect ratio of the dust disk. 



along the curves given by Eq|5]will depend on the surface den- 
sity and temperature T{r) profiles. On the other hand, the 
spatial spread perpendicular to this nominal shape depends only 
on the intrinsic linewidth Av. With sufficient angular resolution, 
we thus expect that Av is only very weakly coupled to any of the 
other disk parameters. The angular resolution should be at least 
high enough to measure the shape described by Eq|5] i.e., mea- 
sure the inclination ;. A more stringent requirement is to resolve 
the transverse size perpendicular to the iso-velocity curves. An 
insufficient resolution would introduce a coupling between Av 
and Z or T as a smaller spatial spread can be compensated for by 
higher temperature (or opacity) to lead to the same flux density. 




Fig. 3. Regions of the disk that yield equal projected velocities 
v b s . The ellipse is the projection of the disk's outer edge. 



3.2. Simple Model 

CS column densities and rotational excitation temperature dis- 
tributions are parameterized using a power law of the radius: 

To derive quantitative information from the data, we used the ra- 
diativ e transfer code adapted to disk DISKFIT (see iPietu et ail / r \ p 
120071 for a detailed description). In a first step (Model A), the Scs(r) 1 £~] 
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T(r) = T \- 



(7) 



As explained by iPietu et all d2007l) . T(r) is the line excitation 
temperature, and it is assumed to apply to all other levels to de- 
rive E(r). We further assume LTE, so that T(r) is now also the 
kinetic temperature. Then, the local linewidth consists of a ther- 
mal component, plus a turbulent contribution 



AV(r) 



i 



2kT(r) 



pm H 



SV tu (r) 2 



(8) 



where \i — 44 is the CS molecular weight and ntn the atomic 
mass unit. We note that AV(r) is a half width at 1/e and that 
the full width at half maximum would be 1.66 times larger. The 
turbulent component is also parameterized by a power law 



(9) 



5V tu (r) = dV I — 



Finally, CS is assumed to be homogeneously distributed in the 
z direction, with a scale height given by the disk midplane tem- 
perature. The density is given by 



«cs = 



■exp(-(z///(r)) 2 ) 



and the scale height is set to 
H{r) = H (r/R h ) L25 



(10) 



(11) 



with H = 16 AU at R h = 100 AU. This definition leads to 
a scale height that is V2 times larger than the H(r) = c JCIk 
convention, where c s is the sound speed and Q.k the Keplerian 
angular velocity. 

All the other geometric disk parameters (systemic velocity, 
inclination and position angle of the disk axis) were taken from 
the an alysis of the CO isotopologues performed by IPietu et all 
(l2007h . We also checked that CS gave consistent results with 
these independent measurements and that the errorbars on these 
geometric parameters did not affect the other derived parameters. 
A summary of the disk characteristics is given in TableQ] 

Synthetic images are generated using the radiative trans- 
fer code DISKFIT, and model visibilities computed for each of 
the observed u, v points. Minimization on the disk free param- 
eters is then made using a x 1 criterion based on the difference 
between the observed and simulated visibilities. This process 
avoids the nonlinearities that result from deconvolution. A mod- 
ified Levenberg-Marquardt method is used to locate the minima, 
with multiple restarts to avoid secondary minima. The linewidths 
we are searching for are not very large compared to the effec- 
tive correlator resolution (0.126 km.s -1 , see Section|2|. Thus, we 
simulated the correlator spectral response by oversampling by a 
factor 4 in velocity space, i.e., using a velocity sampling bin of 
0.02 km.s -1 , followed by a convolution with a kernel mimicking 
the Welch apodization. This kernel extended over 12 consecutive 
channels in the oversampled frequency space. 

Under these assumptions, the CS(3-2) line emission can be 
used to constrain Y,q, p,To,q,dVo and e v . Errobars around the 
minimum are derived from the covariance matrix. A proper 
choice of the pivot values R,,R V ,R S must be made t o mini- 
mize the errors on the power laws (see IPietu et aTll2007l for de- 
tails). With our angular resolution and CS images, R t = R v = 
R s = 300 AU is a good compromise. A degeneracy occurs 
between Eq and Tq if the line is optically thin (for example, 




Fig. 4. Derived nonthermal linewidth dVo at R v = 300 AU as a 
function of assumed kinetic temperature profile. To is the tem- 
perature at R, = 300 AU. The 4 curves correspond to different 
exponents q — (red), 0.2 (green), 0.4 (blue), and 0.6 (cyan). 
Errorbars are ±lcr. 



in the high temperature limit, the emission only depends on 
So/To), but this degeneracy is broken if the surface density 
profile is steep eno ugh to produce an optically thick core (see 
iDutrev et al.l l2QQ7h . Indeed, from a global fit, we find T(r) 
(7.2 ± 0.4)(r/300AUr a63±0 - 09 K. Assuming this temperature to 
be the kinetic temperature, we find an intrinsic nonthermal half 
width at 1/e, 6V,{r) = 0.13 + 0.03(r/300AU)-°- 38±a45 km.s-'. 
The outer radius is found to be 540 ±10 AU. 

The derived excitation temperature is unlikely to differ much 
from the kinetic temperature, because the densities in the disks 
are substantially higher than the CS(3-2) critical densities, even 
at the outer radius. To explain the observed total intrinsic width 

by thermal motions (6V± = ^2kT(r)/pcs>nH, with pes = 44) 
would require a high value of T ^ 50 K at 300 AU. This simple 
analysis thus indicates that pure thermal motions cannot repro- 
duce the observed emission pattern, under the assumptions of 
Model A. 

To strengthen this conclusion, as our main goal is to con- 
strain 6V t (r), we also explored a much wider range of possi- 
ble temperature profiles, varying To and q. Results for dVo as 
a function of assumed Tq are given in Fig|4]for four values of 
q = 0, 0.2, 0.4, and 0.6. In each case, because the exponent e v is 
not well constrained, we assumed e v = q/2, which corresponds 
to a constant Mach number. Figure |4] shows that, unless the ki- 
netic temperature is very high, the nonthermal linewidth must be 
= 0.10 - 0.15 km.s~' to represent the observations. 

The CS(3-2) data set itself suggests low temperatures, but 
the surface density profile is essentially flat (p = 0.14 + 0.20), 
so the constraint may be considered as relatively weak. The tem- 
perature found with C S(3-2) is quite comp atible with the one 
determined from CCH dHenning et al]|2010h and CN and HCN 
(IChapillon et al.ll2012bh . Moreover, the CS(5-4) data, although 
quite noisy, can be used to rule out high temperatures. Figure 
[5] shows the integrated line flux for the CS(5-4) data using the 
best fit parameters found from CS(3-2) as function of To and q. 
Clearly, high temperatures (> 20 K) can be ruled out, since they 
would produce a much brighter CS(5-4) emission than observed. 
This is consistent with the temperature derived from other trac- 
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Fig. 5. Predicted CS J=5-4 line flux from the best fit models 
derived from CS J=3-2 observations, as a function of assumed 
kinetic temperature profile. The 4 curves corresponding to the 
different exponents q = 0, 0.2, 0.4, and 0.6 are essentially de- 
generate. Errorbars are ±l<x. The shaded gray area indicates the 
measured line flux and its +lcr range. Calibration uncertainty is 
not included. 



ers, although subthermal excitation of the CS(5-4) in a warmer 
medium cannot be ruled outQ 

Low temperatures are thus strongly suggested, although not 
firmly proven, by this study. In any case, the uncertainties on 
temperature profile have a negligible influence on the key pa- 
rameters of this study, 5Vq, as well as on pes- However, they can 
affect the derived CS column density by a factor 1 .5-2. 

3.3. A more realistic approach ? 

The assumption of homogeneous vertical mixing (constant abun- 
dance as a function of z) is a strong prior. In practice, chemi- 
cal models predict that molecules are strongly depleted in the 
disk plane, owing to the condensation of volatiles on dust grains. 
Thus, the emission of molecules is expected to come from one 
or two scale heights above the disk plane. As the scale height is 
typically 10 % of the radius, this means the emission comes from 
regions inclined ±Si - 6° from the disk midplane. Because the 
line of sight velocity is proportional to sin(/), these two different 
inclinations result in a broadening of the projected linewidth, 
which is not properly reproduced by the assumption of no ver- 
tical abundance gradient. This is essentially as if the emission 
came from two disks, one inclin ed by ; + Si and the other by 
; - 5i along the line of sight (see lSemenov et ai1l2008l) . To first 
order, the extra broadening due to these inclination differences 
would then be proportional to the rotation velocity, i.e. would 
decrease as y[(r), if Si is constant. The exponent of the nonther- 
mal linewidth law e v (see EqO derived above is indeed quite 
close to the Keplerian exponent 0.5 (but with little significance). 

In practice, Si will vary with radius. On one hand, h(r)/r in- 
creases slowly with r (as as r ). On the other, the height above 

1 This analysis also implicitly assumes a lack of vertical (excitation) 
temperature gradient. On one hand, subthermal excitation is less un- 
likely for the J=5-4 line than for the J=3-2, but on the other, the J=5- 
4 has higher opacities, and can thus probe higher, warmer levels. The 
balance between the two effects cannot be estimated without detailed 
non-LTE radiative transfer models. 



which molecules are found is expected to decrease in the outer 
parts of the disk, because (interstellar or scattered) UV photons 
can penetrate deeper towards the disk plane. This leads to en- 
hanced photodesorption, and brings the molecular layer closer 
to the disk plane. As a consequence, the resulting increase in 
inclination spread is smaller in the outer regions than closer to 
the star. This could also contribute to the apparent decrease in the 
nonthermal broadening found previously by our simple analysis. 

It is thus essential to check whether nonthermal motions are 
also required to explain the observations when CS is confined 
above the disk midplane. Although sophisticated chemical mod- 
els exist, they require a prior knowledge of the H2 density struc- 
ture. This density structure is still very uncertain, and parametric 
chemical models (which would explore a range of possible solu- 
tions) do not yet exist in this respect. 

We thus adopted a second parametrization of the CS emis- 
sion from the disk, Model B, which differs from Model A in two 
ways. The CS surface density is no longer a power law. Instead, 
the disk surface density (dominated by molecular hydrogen) fol- 
lows an exponentially tapered law 

5W(r) = So(^j \xv{-{rlR c f-y), (12) 

which is the self-similar solution of the disk evolution when 
the viscosity is a (time c onstant) power law of radius 
dLvnden-Bell & Pringld[T97l . CS is then assumed to only be 
present when the hydrogen column density from the current 
point towards the disk exterior (perpendicular to the disk mid- 
plane) is lower than some given threshold, S^, the depletion 
threshold. In regions where CS exists, we further assume that 
its abundance is only a power law of radius. The disk is assumed 
to be vertically isothermal. The column density above a given 
height z is then given by 

S K> (r >Z ) = J „(r, M y M =^2|l-erf|^JJ, (13) 

from which the depletion height Zd(f) can be recovered as a func- 
tion of Sd/S(r) 

z d (r)=#(r) erf" 1 (l"f^)- ( 14 ) 
The CS abundance is therefore given by 

Xbs(r)= for z<z d (r) (15) 

Xcs(r) = X° cs (r/R r pcs for z > z d (r). (16) 

In this model, the overall surface density of CS thus follows the 
law 

/ r Y Pcs 

Zcs(r) = X° cs 1 — 1 min(Z d , E H (r)). (17) 

In the modeling, R c and y are fixed paramete rs. We used R c _ 
180 A U and y = 0.55, which were found by iGuilloteau et al.l 
(1201 lh from high-resolution, dual-frequency continuum, and £0 
is also taken from this article to obtain a disk mass of 0.03M o , 
which gives S (100AU) = 9.6 10 23 cm -2 (H 2 molecules). 

Although imposing the density structure could allow non- 
LTE modeling to be performed, considering the large uncertain- 
ties on the real density structural CS is assumed to be thermal- 
ized because densities are in general high enough, even at one to 

2 As stated before, the dust density distribution in disks (and a fortiori 
in the DM Tau disk) is poorly known. For example, the adopted values 
do not consider the variations in the dus t prope rties with radius found by 
IGuilloteau et al.l ( T2011l) . [Andrews et al .1 j2012l> also argue that no simple 
model represents both the gas and dust distributions in TW Hya. 



5 



Guilloteau et al.: Turbulence in disks: the CS tool. 



two scale heights. The remaining free parameters are now X { ^ s , 
Pes, Tq, q, R out , the turbulent width dVo (ey being poorly con- 
strained, we used ey — q/2, a constant Mach number), and the 
depletion threshold 

Although the rationale for the choice of parameters is based 
on behaviors predicted by chemical models, it is important to 
stress that our approach is still purely parametric. In particular, 
we make no a priori link between the thermal structure of the 
disk and the location of molecules. 

The optimal value of was first located through a grid of 
best-fit solutions computed with SV and and To as fixed pa- 
rameters, and X^ s , pes and R out as free parameters. The tem- 
perature was set to the best-fit solution found for the power law; 
similar results were obtained when allowing the temperature law 
to vary. This grid search did not reveal any secondary minimum. 
We then started from the best fit found in this grid search to 
proceed with a global fit of the parameters, including the tem- 
perature law. The best fit results are presented in Table Q] This 
model provides a slightly better fit to the CS(3-2) data than the 
power-law case. 

3.3.1. Turbulent width 

We find a nonthermal component slightly smaller than that de- 
rived in the power-law case, but still significant dVo = 0.11 ± 
0.015 km.s -1 for Ho = 16 AU. A lower value is expected, since 
part of the spatial extent at a given velocity is due to the spread in 
inclinations, which is larger in Model B than in Model A (which 
assumed vertically uniform CS abundances). 

3.3.2. Depletion column density 

The best fit "depletion" column density is = io 2L7±01 cm -2 . 
However, interpretation of this value is not straightforward. 
Changing the scale height from 16 AU (at 100 AU) to 8 or 32 
AU does not significantly affect the derived This indicates 
that l,d is not constrained by the geometry alone. In fact, this 
value is found because the CS(3-2) images are apparently better 
represented by a simple power law for the CS surface density 
than by the more complex radial profiles that can result from the 
hypotheses made in Model B. With ^ 10 22 cm~ 2 and a power 
law for the CS abundance, the CS surface density behaves as a 
power law (of exponent p = 0.4 + 0.2) out to R » 400 AU, close 
enough to the outer radius of the CS distribution. With Ho = 16 
AU, model B also leads to slightly higher (but less constrained) 
temperatures than Model A, around 10 ± 2 K at 300 AU. 

3.3.3. Vertical location 

Although Zrf is mostly constrained by the radial brightness pro- 
file in this case, it does not imply that the vertical distribution 
is not constrained. Starting from the best fit of the radial CS 
distribution, we can actually get an estimate of the system ge- 
ometry by treating the scale height as a free parameter. We find 
Hq = 9+1.5 (see Table [TJ, which is consistent with the tem- 
perature derived from CS, which indicates Hq = 10.4 ± 1.2. The 
required turbulent width increases when the scale height is di- 
minished. The depletion height Zd derived from this best fit is 
plotted in Fig|6l together with the corresponding value of It is 
important to note that the shape of Zd is prescribed by our model 
assumptions. In reality, the distribution of molecules can be sub- 
stantially different, but our analysis nevertheless indicates that a 
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Fig. 6. Top: Total disk surface density E(r) and column density 
Srf(r) above the depletion height Zd(f). Middle: Depletion height 
Zd(r) as a function of radius in the best fit model. The scale height 
H(r) is also indicated (in red). Bottom: as above, but for Zd(r)/r 
and H(r)/r. 



substantial fraction of the CS molecules must reside around one 
scale height. 

3.4. Comparison with chemical models 

A chemical model for sulfur bearing molecules was developed 
bv lDutrev et al.l d201 lb to explain the integrated line profiles ob- 
tained with the IRAM 30-m. The model predicti ons were done 
with the NAUTILUS code dHersant et all l2009t) . which com- 
putes the chemical composition of the gas and interstellar ices 
as a function of time. Bimolecular reactions between gas-phase 
species and interactions with direct UV photons and cosmic-ray 
particles, as well as UV photons induced by cosmic -rays, are 
considered, with a network of 4406 reactions for 460 species. 
Interactions with grain surfaces (including direct and UV in- 
duced dissociations), as well as reactions at the surfaces, through 
1773 reactions, determine the abundance of 195 surface species. 
The NAUTILUS code is run in the vertical ID approximation at 
different radii of the disk (from 50 to 500 AU). More details on 
the model p arameters and the ph ysical structure of the disk can 
be found in iDutrev et al.l (1201 ll) . The observations of CS pre- 
sented injlusy)a])er have b een compared to the Case C for DM 
Tau of lDutrev etaD d201ll see their Tables 4 and 5). 

Although this model does not use the same underlying den- 
sity structure, it is interesting to compare its predictions to the 
results derived from the Plateau de Bure data. Figure [7] shows 
the distribution of CS molecules and the IL density and temper- 
ature vertical profiles for several radii. The range of heights at 
which CS is found corresponds to an IL column density towards 
the outside, Soo(r, z) in the range of 10 22 - 10 23 cirT 2 , except at 50 
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Fig. 7. Predicted distribution of Temperature (top), H2 (middle) 
and CS (bottom) molecule densities as a function of height for 
several radii: 50 (black), 100 (red), 200 (green), 300 (blue), and 
500 AU (magenta). 



AU where the warmer interior leads to enhanced CS abundance 
at larger depths. 

The temperature is not the key parameter defining the 
molecular-rich layer. CS exists at temperatures well below its 
evaporation temperature (30 K), as can be seen in Fig|7] The 
penetration of UV radiation (scattered from the star and coming 
from the ISM) is more important in defining this layer, which is 
limited by photodesorption near the disk plane and photodestruc- 
tion of all molecules near the disk surface. Temperature effects 
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Fig. 8. CS molecule density as a function of (essentially H2) col- 
umn density for several radii: 50 (black), 100 (red), 200 (green), 
300 (blue), and 500 AU (magenta). 



become significant only above 30 K, for which high CS densities 
can be reached because of thermal desorption. 

The enhanced CS surface density within 50 AU from the 
star is not probed by our observations, which are mostly sen- 
sitive to the 100 - 500 AU region. In this region, the predicted 
colum n density increases slightly with radius (see iDutrev et al.l 
1201 ll Fig.4), while the analysis of our observations indicates 
a roughly constant or slightly decreasing column density. The 
chemical model does not include any radial or vertical mix- 
ing, that would be a natural consequence of the observed tur- 
bulent velocities. Vertical mixing will on average decrease the 
altitude at which CS molecules are found, while radial mixing 
will smooth out the CS abundance distribution. Consequently, 
radial mixing would tend to make the slope of the CS surface 
density profile closer to the one of H2 that is expected to vary 
roughly as l/r. This smoothing effect due to mixing can be seen 
in the study of ISemenov & Wiebd d201 ll see their Fig.l 1), who 
investigated the impact of radial and vertical turbulent mixing. 
ISemenov & Wiebd (1201 ll) also find that the CS content can be 
increased by an order of magnitude by mixing. Higher angu- 
lar resolutions (0.1 - 0.3") with ALMA are required to directly 
probe the inner disk region and prove (or disprove) our theoreti- 
cal prediction of high CS abundance here. 

The chemical models predict that CS molecules are located 
between z/r — 0.2 and 0.4 (see Fig|7]i, while our best fit solution 
indicates Zdlr ~ 0.15 - 0.2. When expressed as a function of S^, 
the model predicts the range where CS is abundant reasonably 
well, a few 10 22 to * 10 23 cm 2 as shown by FigJH] The small 
discrepancies suggest that better agreement between the chem- 
ical model and the observations can be obtained by assuming 
larger dust grains and a geometrically thinner (i.e. colder) disk 
in the model. In this case, the mo lecular layer is expected to be- 
come closer to the disk plane (e.g. lAikawa & Nomural2006h . be- 
cause of enhanced UV penetration and lower sticking efficiency. 
Dust settling would also bring the molecular layer in the same 
direction, by minimizing the sticking efficiency above the disk 
plane, but should be less important than grain growth since it 
only concerns the largest grains that offer less surface area per 
unit mass. 
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4. Discussion 



Table 1. Disk parameters and fit results 



The analysis presented above unambiguously indicates that a 
substantial additional line broadening is required besides the 
thermal component and that the derived magnitude of the turbu- 
lent motions is not too critically dependent on the assumptions 
about the location of the CS molecules. 

In the simple power-law approach (Model A), using our 
temperature measurement, a mean molecular weight fi = 
2.30, the derived sound speed is C s (r) - ^kT(r)/ irnin = 
0.16(r/300ALTr ' 3 km.s -1 (using the temperature derived from 
CS for the kinetic temperature), which corresponds to a Mach 
number M - 6V,/C S / V2 around 0.5 and a slightly lower value 
for Model B. 

This result contrasts with the expectations from the a pre- 
scription of the viscosity, where the turbulent width is re- 
lated to <( in a way that depends on the nature of the tur- 
bulence. T he value of 5V ranges between a few aC s to 
~ \[aC s dCuzzi et al. Il200ll). All indirect measurements of 
the a parameter ( Hartmann et al.|[l998t iHueso & Guillodl2005l 
iGuilloteau et al.ll201 II) indicate a values between a few 10 -4 and 
about 0.01, thus predicting much smaller turbulent widths, of at 
most 0.04 km.s 1 if extrapolated near 300 AU. However, these 
measu rements often sample smaller regions. lHartmann et al.l 
(fl998h derive these from accretion rates, and are thus sensiti ve to 
the inner AUs o n ly. Th e results of IHueso & Guillotl (I2005h and 
IGuilloteau et al.l d201 lb are based on viscous spreading of the 
dust disk, and are mostly sensitive to the 100 AUs range. 

iHughes et al.l (I20TT1) obtained a low value for TW Hya 
(< 0.04 km.s -1 ), but up to 0.3 km.s -1 for the Herbig Ae star 
HD 163296 from CO J=3-2 observations. An attractive explana- 
tion of these very different turbulence levels could have been 
the age difference between these two stars, with turbulence 
level decreasing with age as independently suggested by the 
IGuilloteau et all (1201 ll) results. Our DM Tau result does not fit 
well in this scheme, as DM Tau is rather old, too. Differences in 
inclinations between the nearly face-on disk of TW Hya and the 
more inclined ones of DM Tau and HD 163296 could also play a 
role if the turbulence had significant anisotropy. Ho wever, shear- 
ing bo x simulations of MRI induced turbulence bv lSimon et al.l 
( 201 lb show no evi denc e for such anisotropy . On the other hand, 
iFlock et all (l201ll) and iFromang & Nelson! (120061) find small 
anisotropy, the radial velocity fluctuations being about 1.5-2 
times larger than in the azimuthal or vertical directions. 

It is tempting to ascribe the difference between the con- 
straints on a (sampling the disk plane) and our measured tur- 
bulence level (sampling 1-2 scale heights) to the stratification 
of turb ule nce predicted by M HD models. IFromang & Nelsonl 
d2006l) and IFlock et ail (1201 lb show that the velocity dispersion 
increases quickly with height above two scale heights (with the 
c s /Q,k definition, i.e. at y[2H(r) with our definition). However, 
they only reach a Mach number of 0.4 at four scale heights 
(3H(r)), which appears insufficient to explain our result. 

The discrepancy between the apparent linewidths and the es- 
timates of the viscosity may also point towards a weakness in 
our theoretical understanding of the turbulence. The turbulent 
width measures velocity dispersion, while the a viscosity mea- 
sures a transport co efficient . They are related through the effec- 
tive mixing length (lPrandtlHl925h . which may not be a simple 
function of a and H, as is usually assumed in the a model. In 
other words, the aforementioned relation between the turbulent 
Mach number and a may reach its limits in the disk outer parts. 
It is worth pointing out that the TW Hya disk is small (outer ra- 



Geometric 


Adopted 


Fitted 




Parameter 


Value 


Value from CS 




Distance (pc) 


140 






PAD 


65 


65 ±2 




if) 


-35 


-35 ± 1 




Vlsr 


6.08 


6.08 ± 0.02 




V100 (t) 


2.16 


2.17 ±0.10 




M, (M ) 


0.54 


0.54 ± 0.04 




h 


-1.25 






Fitted 


Density Model 




Value 


(A) Power Law 


(B) Tapered Edge 


Note 


x 1 


2468353 


2468336 




H (AU) (a) 


[16] 


9± 1.5 


(1) 


T (K) (b) 


7.2 ±0.4 


8.0 ± 1.3 




q 


0.63 ± 0.09 


0.60 ± 0.20 




Zcs (cm- 2 ) (b) 


5.9 ± 2.5 10 12 




(2) 


X cs (b) 




4.2 ±4.8 10-'° 


(2) 


Pes 


0.13 ±0.20 


0.39 ±0.18 




E rf (citT 2 ) 




a jq21.7±0.1 


(3) 


Rom (AU) 


540 ± 10 


> 580 




dV (km.s" 1 ) (b) 


0.13 ±0.03 


0.12 ±0.025 




ev 


0.38 ± 0.45 


[0.3] 


(1) 



Notes, (t) Rotation velocity (km.s -1 ) at 100 AU, which determines the 
stellar mass M». (a) at 100 AU, (b) at 300 AU. (1) a number between 
brackets [] indicate a fixed parameter. (2) Large errorbar due to strong 
coupling with temperature. (3) Error bar not symmetric; derivation from 
covariance matrix inaccurate. 



dius 230 AU), while those of HD 163296 and DM Tau are large 
(> 600 AU and 900 AU, respectively). The smaller turbulent 
width found in TW Hya is consistent with the idea that viscous 
spreading is the main driver for the disk size. Measurements of 
the nonthermal linewidths for a wider sample of objects with dif- 
ferent disk radii would be important in establishing this potential 
link. 

It is also important to realize that, in all the analysis of the 
molecular line emission, the so-called "turbulent" width is an 
adjustable parameter that is used to fit deviations from the sim- 
ple Keplerian disk model with purely thermal line broadening. 
Departure from these ideal conditions will affect the derived 
"turbulent" width. For example, we have explored in Sect[3]the 
impact of the location of molecules above the disk plane, but 
other effects can exist. Small deviations from Keplerian rotation, 
which are ignored in our analysis, would have no impact on the 
linewidth for a nearly face-on object like TW Hya, but could af- 
fect more inclined objects like DM Tau (i = 35°) or HD163296 
(i = 45°). Our study of the depletion depth in Sectf3]suggests that 
disk warps would, however, need to be quite significant (~ 10°) 
before affecting the required turbulent component. 



5. Summary 

Using the IRAM Plateau de Bure interferometer, we imaged the 
CS (3-2) emission from the disk of DM Tau at ~ 1" resolu- 
tion. Although noisy, CS(5-4) data obtained with the 30-m and 
at 2.5" resolution from the PdBI preclude temperatures exceed- 
ing about 20-30 K. The morphology of the emission indicates 
a significant local linewidth, around 0.14 kms -1 near 300 AU. 
This result appears robust with respect to the assumed location of 
CS molecules above the disk's midplane. While the magnitude 
of this line width is in agreem ent with measurements prev iously 
obtained bv lPietu etal] (120071) and lChapillon et al.ld2012bl) from 
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other molecules like HCO + , CN, HCN or CO isotopologues, the 
relatively large molecular weight of CS (44), combined with the 
temperature estimate and an accurate handling of the spectral 
response of the correlator, leaves no doubt here that this local 
linewidth is essentially due to a nonthermal (turbulent) compo- 
nent. The magnitude of these turbulent motions corresponds to a 
Mach number around 0.3 - 0.5. 

Our data confirm the relatively low temperatures de- 
rived from other molecules such as CN, HCN, and CCH. 
dChapillon et al 11201 2bHHenning et al.ll2010l) . They also suggest 
the CS molecules are located above one scale height from the 
disk plane, in reasonable agreement with predictions from chem- 
ical models. 

This work shows the potential of CS as a turbulence tracer. 
No other mol ecules detected so far ( including the recently re- 
ported HC^N. IChapillon et alj|2012al) is both heavy and abun- 
dant enough to provide a more sensitive tracer. Also, CS offers a 
number of transitions in the mm / submm domain, which probe 
critical densities ranging from a few 10 s to several 10 8 cm -3 . 
This, combined with its expected location in the molecular-rich 
layer, make this molecule a key probe of physical conditions in 
disks. The constraint obtained for the geometric location of CS 
molecules with respect to the disk plane, although still difficult 
to interpret, clearly indicates the potential of ALMA for probing 
different depths. 
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